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abstract 

We investigate the dealer model - an artificial market model based on deter- 
ministic agents both numerically and theoretically. The agents refer to the past 
market prices and changes their ask/bid price. The temporal development of 
the market price fluctuations is calculated numerically. A probability density 
function of the market price changes has power law tails. Autocorrelation co- 
efficient of the changes has an anti-correlation, and autocorrelation coefficient 
of squared changes (volatility correlation function) has a long time correlation. 
A probability density function of intervals between two successive transactions 
follows a geometric distribution. The GARCH type stochastic process is approx- 
imately derived from the market changes of the model in a limit case. We discuss 
two factors of the market price fluctuations and display a relation between the 
volatility of the market prices and a demand-supply curve. We conclude that the 
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power law tails and the long time volatility result from both positive and negative 
heterogenous feedbacks in the agents. 

keywords: artificial market model, GARCH type stochastic process, power-law 
distribution, volatility clustering, econophysics 

1 Introduction 

Economically motivated problems are attracting the attention of physicists, economists 
and engineers. Recently physicists have become interested in these problems 
and established a new field, "econophysics", which has been made tremendous 
progress since 1997. The price fluctuation is one of the most exciting interests in 
econophysic [1, 2, 3, 4]. Statistical properties of market price fluctuations were 
clarified from analyzing financial time series empirically. Specifically we focus 
on two statistical properties of price fluctuations; (1) the probability distribution 
(pdf) of price changes has fat tails [?], and (2) the squared or absolute price 
changes have long time correlation, which is well-known as the volatility clus- 
tering [6]. The volatility clustering is measured by autocorrelation function of 
squared or absolute price changes, the volatility autocorrelation function, which 
has a long tail when the volatility is clustered. 

On one hand stochastic models have been proposed in order to describe these 
phenomena. For example, the truncated Levy flights [5, 7], the random mul- 
tiplicative processes [8, 9, 10] and the GARCH models [11, 12]. On the other 
hand, various agent-based models have been proposed in order to explain the 
basic properties of price fluctuations from the viewpoint of complex dynamical 
systems. These studies are called a agent-based approach [13, 14, 15, 16, 17]. 

The first agent-based model of a market was introduced by one of the au- 
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thors (H.T.) et. al. in order to explain why market price apparently fluctuate 
randomly [18]. It was shown that trading involve a kind of likely nonlinear in- 
teraction among agents in general. A pseudo-random walk fluctuation of market 
price results from the effect of chaos in an artificial market. That is verified in 
the market that includes only 3 agents with very simple deterministic strategies. 
Also, it was pointed out that when the agents have the tendency of following the 
latest market trend, large fluctuations such as crashes or bubbles occur sponta- 
neously. The market model is a variant of the deterministic agent-based model. 

Bak et. al. investigated a stock market model with "rational traders" and 
"noise traders". Noise traders' actions depend on their current volatility in the 
market and imitate ways to buy or sell. Rational traders optimize their own 
utility functions. They emphasize that it is important for noise traders to exist 
in order to have fat tails of the pdf [13]. Lux et. al. investigated a stochas- 
tic multi-agent model with the pool of traders divided into fundamentalists and 
noise traders [15]. They reported that a pdf of returns in their model has fat 
tails. These models assume that a market has a balance of demand and supply 
at a unit time to decide the next market price, and the market price changes de- 
pending on an unbalance of the demand and supply. Unlike the stochastic model 
of Lux stimulative extensions to the basic Minority Game (MG), multi-agent 
market models without noise trader, have been developed and are successively 
investigated [16, 17]. In fact these studies demonstrate the stylized facts (the 
fat tails of the pdf of market price changes and the volatility clustering) and the 
numerical results show that the model contains the GARCH properties. However 
the GARCH process was analytically approximated in these studies. 

Specifically GARCH model are useful tools in econometrics. However there 
are few studies that give the theoretical basis of the GARCH from microscopic 
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viewpoint. The aim of the present article is that we develop an artificial model 
characterized by a few parameters and directly derive the GARCH type stochastic 
process from a price change of the model. In our preceding study we proposed a 
dynamical model of a market having mean field interactions of the dealers [19]. 
We derived a Langevin equation with both positive and negative feedbacks in 
a stochastic manner for market prices and, the pdf of time interval between 
successive trading follows exponential decay. But autocorrelation function of the 
squared change does not have a long time correlation. In the present article we 
introduce heterogeneity of the dealers instead of the mean-field interaction in our 
preceding study. 

The present article is organized as follows. In Sec. 2 we introduce an artificial 
model of a market based on deterministic agents who modify their ask/bid price 
depending on past price changes. In Sec. 3 the statistical properties of artificial 
market price fluctuations from numerical simulations are shown. In Sec. 4 we 
show how the GARCH type stochastic process for the market price changes is 
theoretically derived from the market model in a limit manner. The GARCH type 
stochastic process is a familiar model for time series analysts but its theoretical 
basis in microscopic viewpoint has not been clarified yet. In Sec. 5 we discuss 
the volatility clustering demonstrated from numerical simulation in the model 
comparing with the demand and supply curves of the market. Sec. 6 is devoted 
to concluding remarks. 

2 Market model 

We show a brief explanation of the dealer model - an artificial market model 
with many simple deterministic dealers [20]. The fundamental idea of the model 
is an interaction between agents (dealers) through a common board and histor- 
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ical market prices. We consider a market in which N dealers exchange a single 
commodity, for example, in foreign currency market. We assume that the dealers 
have limited information due to a short time trading. It takes a finite time for 
them to decide to buy or sell. Therefore they have each simple strategy to predict 
the next market price from the limited information. By the same token we do not 
assume that they have utility function since they cannot optimize their action in 
short time. Dealers put a sell/buy order on the market. By a competitive mech- 
anism in the market adequate orders are selected, and the next market price is 
determined from the ordered price. 

As shown in fig. 1 we explain the dealer model for two parts: (1) a market 
rule, which describes how to determine a market price from orders, (2) a strategy 
of the agent, which governs how to order a sell/buy into the market and how to 
modify their ask/bid prices. In the following subsections we explain a market 
rule, a strategy of the agent, initial conditions and parameters. 

2.1 Market rule 

A market has a competitive mechanism for orders. We assume that N dealers 
put their order (sell/buy) to a common board. Here Pi(t) represents an ask/bid 
price of the ith dealer at time t. The ask prices and the bid prices individually 
compete in the market. Namely the maximum buying price and the minimum 
selling price are effective in the market. Thus the condition for a trading to occur 
is given by the following inequality, 

max {p(t)} > min {p(t)}, (1) 

tor all buyers for all sellers 

where the right hand side represents the maximum bid price for all the buyers, 
and the left hand side represents the minimum ask price for all the sellers. As- 
sume that the market price P(t) is determined as an arithmetic mean of the buy 
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price and the sell price when trading occurs. Otherwise the last market price is 

maintained. Namely, 
1 

P(t) 



(max{p(t)} + min{p(i)}) (max{p(t)} > min{p(t)}) 

^ for all buyers for all sellers for all buyers for all sellers 



P(t - 1) (max{p(t)} < min{p(t)}) 

for all buyers for all sellers 



(2) 



2.2 Strategy of the dealer 

In general a dealer determines his/her action from several causes. The causes are 
separated into "endogenous" factors and "exogenous" ones. The endogenous ones 
mean what has happened in the market and contains the historical market prices 
and a rumor in the market. The exogenous ones mean what has happened out- 
side market, a balance of domestic commerce or international commerce among 
countries and so on. Here we consider that the dealer has to make a decision 
based on the market price fluctuations since they regards only prices for a short 
period and determine their action and ask/bid price. 

Here we consider a trade between a seller and a buyer. Suppose that the both 
must make a trade. Then the seller decreases his/her sell price a little bit unless 
he/she can make a trade. On the other hand the buyer increases his/her buy 
price a little bit. Repeating this they find a satisfying exchange price. 

The same is consistent in the market. Sellers go on decreasing their sell 
price a little bit until they can sell and buyers go on increasing their buy price 
until they can buy. From the assumption a modification of the ask/bid price is 
negative/positive for sellers/buyers. The temporal development of the ask/bid 
price can be described by, 

p i (t + l)=p i (t) + a i (t)D i (t), (3) 

ai(t) represents the i'th dealer's modification per a unit time step, and Di(t) = —1 
when the ith dealer is a seller, and Di(t) = 1 a buyer. 
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Moreover we assume that a value of modification depends on the past market 
prices due to limited information. The assumption means that the next action 
of dealers is only affected by the historical data of market prices. 

a l (t)=a(P(t),P(t-l),...) (4) 

We assume that the dealers are sensitive to price changes rather than exact prices. 
One of the simplest modification algorithms is a linear inner product of dealer- 
dependent coefficients and past price changes. Namely the modification a^t) is 
described by 

T 

Oi(t) = |1 + Ci,s< ^Pprev{.s')\a U (5) 

s'=l 

where AP pre „(s') denotes the s'th change of the past market price. It is obvious 
that AP prev (l) is the latest market price change. q jS ' represents coefficients of 
the ith dealer for the last s'th market price change, a« is a positive coefficient. 

For simplicity we consider the case of T = 1. Then a rule to modify dealers' 
expectation price is written by 

Pi (t + 1) = Pi (t) + |1 + c i AP prcw (l)|a i A(0, ( 6 ) 

where q is a coefficient, which corresponds to accelerating and deaccelerating 
his/her modification of expectation price depending on the latest price. After the 
large price change each dealer goes up or goes down quickly his/her expectation 
price in order to make his/her portfolio balanced as soon as possible. 

We assume that the dealers open a sell/buy position till they exchange. They 
determine whether they open the sell/buy position after they have made a trade. 
We assume that the sellers want to go on selling when the market price goes up 
and that the buyers want to go on buying when the market price goes down. 

Therefore after trading the commodity the two dealers who have exchange 
open either a sell position or a buy position depending on the past price change 
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AP prev . If the market price goes up then the dealers may expect to make a profit 
from selling the commodity. Hence we assume that they open the buy position, 
i.e., Di{t + 1) = —1, when AP prev > 0. In striking contrary we assume that they 
open the buy position, i.e., D^t + 1) = 1 when AP prev < in order to make 
a profit from future selling the commodity. Furthermore they open either a sell 
position or a buy position at the same probability p — 1/2 when AP prev = 0. 

Moreover we assume that dealers have difference between the market price 
and the retried ask/bid price, which denotes Aj(= P(t) — Pi(t)). In a mind of a 
dealer it is clear that the retried ask price is greater than the market price. On 
the other hand the retried bid price is less than the market price. Therefore the 
ask/bid price of the ith dealer at t is written by, 

Pi (t)=P(t)-D(t + l)A i (7) 

For simplicity Aj = A for the all dealers. 

2.3 Initial conditions and parameters 

Each dealer has two coefficients: Oj and q. and q are initially given by uniform 
random numbers distributed in interval [0, a*] and [— c*,c*], respectively. These 
coefficients exhibit dealer's personality and are fixed throughout a numerical sim- 
ulation. Each dealer has two variables: Di(t) and Pi(t). Initial condition -D«(0) 
is given by either +1 or —1 randomly, i.e. a probability for Dj(0) to be 1 is 1/2, 
and a probability for -Dj(O) to be —1 is 1/2. Then Pi(Q) is given by P — A when 
Di(0) = 1 and Po + A when D{(0) = —1, where Pq is an initial market price fixed 
as P = 120.0 throughout simulations. 

We assume that the latest price change AP prev is initially zero. The dealer's 
rule is deterministic except initial conditions. This model has four parameters; 
amount of the dealers N, a* for a,, difference between the ask/bid price and the 
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market price A and c* for q. Finally we summarize the model parameters in Tab. 
1. 

3 Numerical simulation 

Fixing parameters N = 100, A = 1.0 and changing c* and a*, we numerically 
simulate the dealer model introduced in the above section. Figs. 2 and 3 show a 
typical example of time series of market prices P(t) and of price changes AP(t) = 
P(t)—P(t—1). The market prices and their changes apparently fluctuate although 
the model is completely deterministic. The reason is high-dimensionality of the 
system as the freedom of the system is proportional to the agent-number. 

Let r s denote the time when the sth trade occurs. Fig. 4 exhibits a conceptual 
illustration of market price fluctuations. The time series may be characterized by 
a price change and a time difference between successive trading. P(t s ) represents 
a market price at t s . Ap s denotes a price change at r s , namely, Ap s = P(t s ) — 
P(t s -i), an d n s a time difference between successive trading, n s = r s — r s _i. 

Fig. 5 displays semi-log plots of pdfs of Ap s for various values of a* and c*. 
The pdfs seem to have fat tails and their tails depend on the values of both a* 
and c*. The corresponding cumulative distribution function (cdf) is defined as, 



where f(x) is a pdf. If the pdf follows power-law distribution the corresponding 
cdf is given by 



where (3 is a power law exponent (f3 > 0), which is estimated from a slope of the 
cdf in the double-log scale. As shown in fig. 6 the power law exponent depends on 
the value of both c* and a*. For large a* and c* the power law exponent becomes 





F(> \x\) oc \x 



(9) 
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small. Moreover let be k = a*c*. Then the power law exponent (3 is a function 
of k as shown in fig. 7. We discuss why the power law exponent depends on k in 
the next section. 

Typical example of time series of n s and its pdf are shown in fig. 8, respec- 
tively. The pdf can be approximated by a geometric distribution, 

Win) =p(l-p) n (10) 

This means that occurrence of trading is fully random. By numerical fitting in fig. 
8 (right) we obtain p = 0.36. From analysis of high-frequency financial data it is 
clarified that the pdf trading intervals at the same time is fitted by exponential 
distribution [21]. 

4 Stochastic approximation as the G ARCH pro- 
cess 

Here we consider why the pdf of the price changes has fat tails and that its tail 
index depends on the parameter both a* and c*. We will show that the market 
price is approximately dominated by the GARCH process through this section. 

As shown by the numerical simulations the market price fluctuates discontin- 
uously. Here we consider Ap s , a change of the market price on the sth transaction 
as shown in fig. 4. Let M s denote the buying price at the sth transaction, and 
m s the selling price as shown in fig. 9. From the definition P(t s ) = \{M S + m s ) 
and Ap s = P(t s ) - P(r s _i), we get 

A Ps = ^(M s + m s ) - ^(M s _! + m s _i) 

= -(M a -M i _i) + -(m s -m a _i). (11) 

From eq. (6) the next buyer at r s _i adds |1 + c J Ap s _ 1 |aj into his/her bid price 
n s times until he/she can trade, and the next seller also subtracts |1 + CiAp s ^i\a,i 
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from his/her ask price n s times. Hence the first term and the second term in eq. 
(11) are given by 

M s -M s _x = \1 + CjAps^lajUs-! -K 9 , (12) 
m s -m s _i = fc s - |1 + c i Ap s _i|a i n s _i, (13) 

where the subscript j represents the dealer who gives the highest bid price at the 
time t s , and i represents the dealer who gives the lowest ask one. K s represents 
a difference of the bid prices between the (s — l)th buyer and the sth buyer at 
r s _i, and k s of the ask prices between the (s — l)th seller and the sth seller. 
Substituting eqs. (13) and (12) into eq. (11) yields 

Ap s = i|l + CjAps-^ajTis-! - i|l + CiAps-^aiUs-! + ^(k s - K s ). (14) 

Since the dealer number N is large we assume that {c} and {a} are mutu- 
ally independent stochastic variable, which is uniformly distributed in intervals 
[-c*,c*] and [0,a*]. Therefore we have (c) = 0, (c 2 ) = c* 2 /3, (a) = a*/2 and 
(a 2 ) = a* 2 /3. Moreover since the sellers and the buyers are symmetric k s — K s is 
also symmetric. Actually it is easily confirmed by a numerical simulation. Fig. 
11 (left) is a typical example of time series of k s — K s . It is obvious that k s — K s 
is symmetric. Therefore (k s — K s ) = 0. 

By taking conditionally averaging eq. (14) over dealers' indices i and j we 
have, 

(A Ps ) = 0. (15) 

By taking square of eq. (14) and averaging over dealers' indices i and j under 
the condition that Ap s „i is realized, we get 

(Ap 2 ) = \{{k s - K s f) + + cA Ps ^){a 2 ){nl) 
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+ ^(|l + cA^ 1 | 2 )(a 2 )(^)-i(|l + cA^ 1 |)(|l + cAp s _ 1 |)(a) 2 (^) 
- 7<(*. - K s ) 2 ) + + cA Ps _ 1 | 2 )(a 2 )(n 2 ) - \(\\ + cA Ps _ 1 |) 2 (a) 2 (n 2 ) (16) 



The first term in the left hand side of eq. (16) is not constant. k s — K s has 
large fluctuations. Fig. 11 (right) is the pdf of k s — K s . The pdf is the same 
distribution as Ap. Hence we assumed that ((k s — K s ) 2 ) is proportional to a 
conditional variance of Ap s _i, namely 

\((k s -K s ) 2 )= v (A P 2 _ 1 ), (17) 

where rj is a positive coefficient. The second term in the left hand side of eq. (16) 
is calculated as, 

+ cA Ps _ 1 | 2 )(a 2 )(n 2 ) = ^(n 2 ) + ApU(n 2 s ) . (18) 

The third term in the left hand side of eq. (16) is calculated as, 

i(|i+cA P ,_ 1 |) 2 (o) 2 ( n ;> = 



2 { %(n 2 ) (|Ap s _i| < 1/c*) 

(19) 



Therefore by using a s = (Ap 2 ) — (Ap s ) (16) is described as, 

2 ( U! + rjaU + Pl ApU + (\Ap s ^\ > 1/c*) 

^ I "2 + WU + P 2 Ap 2 _i (lAp^l < 1/c*) ' 1 J 

where u x = ^a* 2 (n 2 }, Pl = ^ 2 (n 2 s }, 7 = ^(n 2 ), w 2 = ^a* 2 (n 2 ) and p 2 = 
±n 2 (n 2 ). From eq. (10) we obtain (n 2 ) = {2 ' p) ^~ p) « 2.92 

For large |Ap s _i| the term 7-^2 — of eq. (20) is almost zero. Then eq. (20) is 
approximated as a GARCH(1,1) process. The probability density function of the 
GARCH(1,1) process has fat tail [22, 23]. The power law exponent is a function 
of <Ti and i]. Namely the power law exponent depends on k. Furthermore if k is 
small then the second term of eq. (20) vanishes, and we obtain a 2 = 0J2 + Ws-i- 
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The result of the iteration is given by, 

„J _ + (,_«,). (21) 

Namely the variance of Ap s is nearly constant. Then p s behaves similarly to a 
Brownian motion. 



5 Discussion 

5.1 Correlation function 

We calculate an autocorrelation coefficient of the market price changes and one of 
squared changes. The volatility clustering is the well-known fact that the squared 
changes of the market price are clustering in financial data. It is indicated that 
the volatility clustering is related to a long time correlation of the market [6]. 
These two autocorrelation coefficients are defined by 

R (i) _ (Apt+s^Pt) - (Ap t+s )(A Pt ) 

s ~ <A^>-(A Pt > 2 ' {22) 

p( 2) (Ap 2 + ,Ap 2 )-(Ap 2 +s )(Ap 2 ) , , 

3 " (Arf> - (Ap?>2 • [ 6) 

As shown in fig. 10 (left) the autocorrelation coefficient is negative at small s. It 
means that the price change tends to move to an opposite direction of the last 
price change. On the other hand autocorrelation coefficient of the squared change 
in fig. 10 (right) has a long time correlation. It is interesting that the volatility 
of price changes has long time correlation although each dealer just depends on 
the latest price change. Actually it is reported that the volatility ocorrelation 
function of the GARCH(1,1) process decreases exponentially [24]. However that 
is only true when a ARCH(l) parameter plus a GARCH(l) parameter is less than 
or equal to 1. When the ARCH(l) parameter plus the GARCH(l) parameter is 
greater than 1 it is impossible to derive the autocorrelation function in the same 
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analytical way. From numerical simulations of the pure GARCH(1,1) process 
we confirm that the autocorrelation function has long tail when the ARCH(l) 
parameter plus the GARCH(l) parameter is greater than 1. Specifically its long 
tail is significant when the GARCH(l) parameter is near 1. Hence we can explain 
from this property that the market price fluctuations of the dealer model exhibits 
the volatility clustering when r\ + pi > 1 and rj is near 1. 

5.2 Demand and Supply curve of the model 

In the model introduced in the article we treat the case that demand and supply 
are automatically balancing, so that the total number of both sellers and buyers 
is conserved. We introduce time- dependent density of buyers and sellers for a 
market price P at time t, d(P,t) and s(P,t), respectively [25]. The cumulative 
frequency distributions of sellers and buyers D(P,t) and S(P,t) are respectively 
defined by 



Of course in practical markets numbers of sellers and buyers can be observed 
partially. However in the numerical simulation one can calculate them from all 
the dealers' variable. Here we consider relation between volatility and a demand- 
supply curve. Fig. 12 shows D(P,t) — S(P,t) in high volatility regime (A) and 
in low volatility one (B). It is found that a slope of D(P,t) — S(P,t) in the high 
volatility regime is gradual. On the other hand one in the low is rapid. This 
means that a slope of D(P,t) — S(P,t) is related to the volatility. Hence the 
volatility clustering is attributed to gradually changing of the demand-supply 
curve. 




(24) 



roc 

S(P,t) = / s(P',t)dP' 



Jp 



(25) 
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5.3 Amount of sellers and amount of buyers 

In the model the seller and buyer who have exchanged commodity open his/her 
position depending on the latest price change after they make a trade. We assume 
that the sellers go on selling when the market price goes down and that the buyers 
go on buying when the market price goes up. Let be Ns and Nb represent an 
amount of the sellers and one of buyers. Then N s — N B fluctuates around zero. 
The reason is as follows: When N$ — Nb > the market price tends to go down. 
When Ns — Nb < it goes up. However when the market goes down Ns decreases 
and N B increases. N s — N B becomes less than zero. When the market price goes 
up Ns increases and Nb decreases. Ns — Nb becomes greater than zero. If the 
sellers/buyers go on selling/buying when the market price goes up/down then a 
crash/bubble occurs. Throughout both the GARCH process and the demand and 
supply we think that there are two reasons why the market price fluctuates. One 
results from rules of a dealers' prediction depending on the past market price. 
Another results from unbalancing of demand and supply. 

6 Conclusion 

We introduced the dealer model in which a commodity is exchanged such as 
a foreign exchange market. The model generates the time series statistically 
similar to real financial one although the model has a few parameters. This 
means that the model contains primitive factor for the market price fluctuations. 
The probability density function of the artificial price changes has power law 
tails. The autocorrelation function of the changes has anti-correlation for a few 
ticks. This implies that the market price changes tend to move in opposite. A 
squared autocorrelation function (volatility correlation function) has a long tail 
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although the dealers depend on the latest price change. We derived the GARCH 
process for price changes under the assumption that the dealer homogeneously 
trades irrespective of their coefficients i.e. we can take average of Ap"j, over all the 
dealers' indices. From the GARCH(1,1) limit we explained the power-law of the 
probability density function and the long time correlation of volatility correlation 
function. The long time correlation of volatility can be also seen in a demand- 
supply curve. The slope of the demand-supply curve is gradual in high volatility 
regime. In contrary the slope of demand-supply curve is rapid in a low volatility 
case. We guess that the market price fluctuations result from both the dealers' 
prediction depending on the past market prices and unbalancing of demand and 
supply. 

The relation between microscopic behavior of dealers and macroscopic feature 
of a market will be bridged by statistical mechanics. Then the intriguing prob- 
lem is to clarify what kind of microscopic interaction between dealers generate 
statistical feature of macroscopic variables. 

Acknowledgements 

One of the authors (A.-H. Sato) thanks T. Munakata and M. Aoki for stimulative 
discussions and useful comments. 

References 

[1] H.E. Stanley L.A.N. Amaral, P. Gopikrishnan and V. Plerou, Physica A, 
283 (2000) 31-41. 

[2] M. Ausloos, Physica A, 285 (2000) 48-65. 



16 



[3] X. Gabaix, P. Gopikrishnan, V. Plerou and H.E. Stanley, Nature, 423 (2003) 
267-270. 

[4] X. Gabaix, P. Gopikrishnan, V. Plerou and H.E. Stanley, Physica A, 324 
(2003) 1-5. 

[5] R.N. Mantegna and H.E. Stanley, An Introduction to Econophysics: Corre- 
lations and Complexity in Finance, Cambridge University Press, Cambridge 
(1999). 

[6] M. Pasquini and M. Serva, Physica A, 269 (1999) 140-147. 
[7] H. Nakao, Physics Letters A, 266 282-289. 

[8] H. Takayasu, A.-H. Sato and M. Takayasu, Physical Review Letters, 79 
(1997) 966-969. 

[9] D. Sornette, Physica A, 250 (1998) 295-314. 

[10] A.-H. Sato, H. Takayasu and Y. Sawada, 61 (2000) 1081-1087. 

[11] R.F. Engle, Econometrica, 50 (1982) 987-1007. 

[12] T. Bollerslev, R.Y. Chou and K.F. Kroner, Journal of Econometrics, 52 
(1992) 5-59. 

[13] P. Bak, M. Paczuski and M. Shubik, Physica A, 246 (1997) 430-453. 

[14] N.F. Johnson, S. Jarvis, R. Jonson, P. Cheung, Y.R. Kwong and P.M. Hui, 
Physica A, 258 (1998) 230-236. 

[15] T. Lux and M. Marchesi, Nature, 397 (1999) 498-500. 

[16] D. Challet, M. Marsili and Yi-Cheng Zhang, Physica A, 294 (2001) 514-524. 

17 



[17] P. Jefferies, M.L. Hart, P.M. Hui and N.F. Johnson, The European Physical 
Journal B, 20 (2001) 493-501. 

[18] H. Takayasu, H. Miura, T. Hirabayashi and K. Hamada, Physica A, 184 
(1992) 127-134. 

[19] A.-H. Sato and H. Takayasu, Physica A,250 (1998) 231-252. 

[20] A.-H. Sato and H. Takayasu, Empirical science of financial fluctuations - 
The Advent of Econophysics (Ed. by H. Takayasu), Springer- Verlag, Tokyo, 
(2002) 171-178. 

[21] M. Takayasu and H. Takayasu, Physica A, 324 (2003) 101-107. 

[22] D. Ghose and K.F. Kroner, Journal of Empirical Finance, 2 (1995) 225-251. 

[23] P.A. Groenendijk, A. Lucas and C.G. de Vries, Journal of Empirical Finance, 
2 (1995)253-264. 

[24] Z. Ding and C.W.J. Granger, Journal of Econometrics, 73 (1996) 185-215. 
[25] H. Takayasu and M. Takayasu, Physica A, 269 (1999) 24-29. 



Table 1: Model parameters; N, a*, A and c*. 



parameter 


detail 


N 


amount of the dealers 


a* 


a prediction coefficient Oj 


A 


difference between retried ask/bid price and the market price. 


c* 


for a prediction coefficient q 
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market 



market price 
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buy 
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Figure 1: The conceptual illustration of the dealer model. The inputs of the 
market are orders from dealers. The output of the market is a market price. The 
input of an agent is a sequence of past price changes. The output of an agent is 
a sell price or a buy price. 
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Figure 2: A typical example of time series of market prices at N — 100, a* = 0.01 
and A = 1.0 for c* = O.O(right) and 375.0(left). 
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Figure 3: A typical example of time series of price changes at N = 100, a* = 0.01 
and A = 1.0 for c* = O.O(left) and 375.0(right). 
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Figure 4: A conceptual illustration of a market price fluctuation. A trading 
occurs at random. t s represents the time when the sth trading occurs. P(t s ) 
exhibits a market price at t s , Ap s a market price change at r s , namely Ap s = 
P(t s ) — P(r s _i), and n s a difference between r s _i and r s , n s = r s — r s _i. 
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Figure 5: Semi-log plots of the probability density functions of the price changes 
Ap s . We fix N = 100, A = 1.0 and a* = 0.01, and change a* (left). We fix 
N = 100, A = 1.0 and c* = 100.0, and change a* (right). 
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Figure 6: Log-log plots of the cumulative distribution function of the price 
changes Ap s . We fix = 100, A = 1.0 and a* = 0.01, and change a* (left). 
We fix N = 100, A = 1.0 and c* = 100.0, and change a* (right). 
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Figure 7: Log-log plots of the cumulative distribution functions of the price 
changes Ap s at k — 3.8 and k = 3.7. It shows the cumulative distribution 
functions at three different parameter sets for each k. 
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Figure 8: A typical example of time series of an interval between successive 
trading n s (left) and its probability distribution function in semi-log scale (right). 
Parameters are fixed at iV = 100, a* = 0.01, A = 1.0 and c* = 375.0. 
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Figure 9: A conceptual illustration for explaining dealers interaction between 
t s _i and t s . M s represents a buy price at t s , and m s a sell price. The Next buyer 
at r s _! adds |1 + CjAp s _i\ into his/her bid price n s times, and the next seller 
subtracts |1 + CjAp s _!| from his/her bid price n s times. 
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Figure 10: Autocorrelation coefficient of price changes (left) and Autocorrelation 
coefficient of squared price changes (right) are calculated from temporal develop- 
ment of price changes. Parameters are fixed at N — 100, a* = 0.01, A = 1.0 and 
c* = 375.0. 




Figure 11: A typical example of time series of k s — K s at N = 100, a* = 0.01, 
A = 1.0 and c* = 375.0 (left). The pdf of k s — K s and Ap s at the same parameters 
(right). 
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Figure 12: A difference between cumulative frequency distribution D(P,t) and 
S(P,t). Capital A exhibits a high volatility regime, and B a low volatility one. 
Parameters are fixed at A = 100, a* = 0.01, A = 1.0 and c* = 375.0. 
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